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Quantum devices based on Josephson effect in superconductors are usually described by a Hamil¬ 
tonian obtained by commonly used canonical quantization. However, this recipe has not been yet 
rigorously justified. We show that this approach is indeed correct in certain range of parameters. We 
find the condition of the validity of such quantization and the lowest corrections to the Josephson 
energy. 


The Josephson effect [I; is one of most fundamental 
phenomena occurring in physical systems with macro¬ 
scopic quantum effects. When biased by external capaci¬ 
tance, inductance or current source, it becomes the basis 
for the description of many superconducting quantum 
devices such as Cooper pair boxes, flux and phase super¬ 
conducting qubits M- It has been demonstrated that 
such structures can be promising elements for quantum 
computing [4j and violation of local realism (Bell test, 
although only local at present) [5j. The detailed theoret¬ 
ical description of such a system becomes more and more 
necessary for its highly complex and demanding experi¬ 
mental applications. 

The simplest, and to date most popular and phe¬ 
nomenological description of the quantum devices based 
on Josephson effect is derived from current-phase rela¬ 
tion, phase-voltage relation. The quantum properties of 
Josephson junction are extracted by using conjugated ob¬ 
servables for example between number particles (charge) 
and phase 0. The energy stored in capacitor decides 
whether we deal with Cooper pair box, flux or phase 
qubit. On the other hand, the whole system is already 
originally quantum, with superconductor (conventional, 
s-wave) described at least by standard Bardeen-Cooper- 
Schrieffer (BCS) theory |J[ [8] so there should be no need 
to quantize it again. The canonical quantization pro¬ 
cedure is just a shortcut between the full quantum de¬ 
scription and the approximated model [5]. Certainly one 
should be able to justify this procedure and find the range 
of its validity. A partial justification has been developed 
30 years ago [TO] , where the effective action of the Joseph¬ 
son junction has been found and - in path integral pic¬ 
ture of quantum mechanics - the canonically quantized 
Hamiltonian is revealed in the tunneling limit. However, 
that approach does not give an answer to a very rele¬ 
vant question of the range of validity, except qualitative 
statement that phase fluctuations should be smaller than 
the superconductor energy gap. The point is that tun¬ 
neling limit (low transmission through the junctions) is 
not a sufficient condition as both transmission, number 
of modes and capacitance play the role. If one simply 
takes tunneling limit at constant capacitance then also 
whole useful quantum properties are lost. Phase fluctu¬ 
ations are indeed intuitively of the order of capacitance 


energy but one has to settle it quantitatively. 

In this work, we derive the approximate Hamiltonian 
(corresponding to the one obtained by canonical quanti¬ 
zation) straight from BCS Hamiltonian, using standard 
perturbative approach with both transmission and ca¬ 
pacitance treated as small parameters. By introducing 
phase-dependent eigenstates, the formalism takes much 
analogy from adiabatic and nonadiabatic transitions. We 
find that in the lowest order, which can be called adi¬ 
abatic approximation CD, one indeed finds the desired 
simple effective Hamiltonian. The next order modifies 
Josephson energy and ratio of this modification gives 
the range of applicability of the adiabatic approxima¬ 
tion. Interestingly, the modification contains an infinite 
term, which is fortunately phase-independent and so it 
has no physical relevance, and appears also in the non¬ 
superconducting case. The phase-dependent part can be 
calculated analytically also for large transmission, but it 
is rather technical so the details are presented in Supple¬ 
mental Material [12] . Higher order terms become rather 
complicated, but do not contain anymore infinity. 

Let us begin with the sketch of the common canonical 
quantization of an electric circuit with Josephson effect 
[B]. The system depicted in Fig. [T] consists of the Joseph¬ 
son junction (JJ), effective capacitance C and bias induc¬ 
tance L or current Iq. The total Hamiltonian will be the 
sum of the contributions of all these elements. We have 
Hq = Q 2 /2C = CV 2 /2, where Q is the charge and V is 
the voltage between 1 and 2. The voltage can be repre¬ 
sented by phase (j) change ratio V = hef/e (2cj> = 3>i — $ 2 
for the respective superconducting phases $ 1 , 2 )- Anal¬ 
ogously, Hl = LI 2 /2 = (h(j)/e) 2 /2L with the current 
I = Q. Alternatively, one can consider current bias / 0 
giving the term Hi = Ighcf/e (Jo is then the current 
between 1 and 2 in absence of other elements). The clas¬ 
sical picture leads to relation Q = I = edH/hdcf. On 
the other hand the Josephson current is calculated as 
Ij = Ii sin(2<^) PQ. Using the above classical equation, 
this leads to Josephson energy Hj = — ehl\ cos(2</>)/2. 
Now, the canonical quantization means that we take 
the Hamiltonian H = He + H^/i + Hj as a function 
of conjugate variables <j> and Q and replace everything 
by their quantum counterparts, including conversion of 
Poisson bracket into commutator between <j> and Q. It 
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FIG. 1: (a)The effective electric circuit of the Josephson 
junction (JJ) biased by capacitance C, inductance L and/or 
current Jo with distinguished parts 1,2. (b) The energy 

Hl/i + Hj with the quantum operation range in the inset, 
(c) Quantum energy levels in the potential landscape 


essentially means Q ed/id<j). We stress that the 
canonical quantization of the LIC part {He + H^/i) 
is already well justified by quadratic form of the elec¬ 
tromagnetic action [T3], whose Lagrangian density has 
the form £ = eo\E\ 2 /2 — \B\ 2 /2p, 0 — pV + A ■ j, where 
E = —\7V — d t A,B = V x A,V,A,p,j,eo,po denote electric 
and magnetic field, scalar and vector potential, charge 
density and current, electric and magnetic permittivity, 
respectively, and linear relations between V and E, Q 
and p 1 I and j. Certainly, one has to assume ideal mod¬ 
els of capacitance and inductance, which will surely break 
down at high frequencies or voltages (e.g. leakage cur¬ 
rents). Nevertheless, the useful range of quantum oper¬ 
ation parameters is well within the regime of validity of 
the assumed model of the LIC part. One usually consid¬ 
ers the classical turning points defined by E ~ Hx + Hj , 
where Hx can be approximated by e 2 cj)Q(j) / Le 2 with the 
exact crossing point <j>o, equivalent to current bias term 
Hj , see Fig. Therefore the most suspicious part 

of this procedure is not the quantization of LI or C 
terms but JJ. On one hand, the Josephson effect is cru¬ 
cial for the appearance of quantum levels in the potential 
landscape, see Fig S'- On the other hand, the effect is 
no longer quadratic and originally derived straight from 
BCS theory plus tunneling. It is intuitive that going 
back-and-forth between quantum and classical descrip¬ 


tion is not always accurate and may lead to losing some 
corrections. To find this accuracy is the main goal of this 
work. 

We shall model the full Hamiltonian of the supercon¬ 
ductor (s-wave) by the standard bulk BCS part written 
in the Bogoliubov-de Gennes framework 0 Hbcs j plus 
tunneling term Hx- We shall use the relation between 
second and first quantized Hamiltonian H = F with 
anticommutation {’F'1,’1'} = I and {\F,\F} = 0 (here \F 
can be multi-component). The first quantization Hamil¬ 
tonian h = Hbcs + hr is defined in a one-dimensional 
bispinor basis x, (1,2) x {e,h), where x is the (real) po¬ 
sition in the superconductor 1 or 2 and e/h denotes elec¬ 
tron/hole. In the standard BCS approximation (small 
superconducting gap), for a single channel, we have 


h-BCS 


f-id x 

\ A 


A 

id. 


, h eT = -tS{x ) 


eh 




(i) 

and hhT = —h e x{4> —> —<j>) where hscs is the same in 
subspaces 1 and 2 while hx has different phase in electron 
and hole subspace. Here A is the superconducting gap 
and position is represented in time units divided by K. 
Normal scattering matrix is then 






( 2 ) 


with t = sin(r/7i) and r = yl — t 2 ,Sh = Sj . One can 
add the overall phase e la ! h to the scattering matrix by 
modifying hx h? + a5(x ) (accounting interface prop¬ 
erties) but it will not change any of our results and hence 
we can safely disregard it. We stress that putting (j) into 
the tunneling term instead of varying phases at A in the 
BCS part is only matter of the gauge choice. In our gauge 
the phase 4> corresponds simply to the vector potential 
across the junction J) eA ■ df/H so that superconducting 
phases remain constant all the time. 

Now one can calculate the ground energy (we focus on 
zero temperature case) of the whole Hamiltonian assum¬ 
ing Fermi level at 0 (all states with negative energy e are 
occupied). The Josephson effect follows from the fact 
that this energy is phase-dependent. The spectrum of 
^bcs consists of continuum positive and negative parts 
separated by the gap |e| > A Introducing hx mod¬ 
ifies the wavefunction of these (now scattering) levels 
but not their energies. However, hx allows for new, so- 
called Andreev bound states (ABS) [13] , with energies 
inside the gap |e| < A and hence localized at the junc¬ 
tion point (x = 0). These energies (see detailed calcula¬ 
tion in Supplemental Material |12| which is anyway quite 
straightforward m) are equal e± = ± A \fl — t 2 sin 2 _d> 
and, for small t, ~ ±A At 2 sin 2 </>/2, see Fig. [2] 
Taking Hj — e_, and summing over independent chan¬ 
nels, Ylti as usually JJ contains many, we have the de¬ 
sired canonical quantized scheme. Note also that the 
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FIG. 2: Energy spectrum of JJ, with the continuum parts 
above +A and below —A and two phase-dependent ABSs in 
the gap. Possible excitations of the ground state: (A) be¬ 
tween ABSs, (B) between an ABS and continuum, (C) be¬ 
tween lower and upper continuum. 

Joseplison current can be justified as I = dedHj/hdcj) by 
means of full counting statistics pbj, where (e lx< ^^) = 

^ e itH(4,-hxl2e)lh e -itH{ < t > +hx/2e)/h^ j n the ground sta te, 

which follows from Keldysh formalism and is in fact 
equivalent to Josephson’s original derivation [I]. 

The above presented analysis, although leads to the 
desired result, is incomplete. It is not only the en¬ 
ergy that changes with the phase 4> but also the whole 
wavefunction. One has to ensure that correction from 
these changes is sufficiently small in the quantum opera¬ 
tion regime, defined by e 2 /C ~ The argument 

that wavefunction changes become small for t —> 0 (i.e. 
•C e 2 /A C) is insufficient because all relevant quantum 
properties (energy levels) disappear in this limit. To find 
rigorously the range of validity of the canonical approx¬ 
imation we have to calculate the effect of wavefunction 
changes. 

We shall treat the problem by standard perturbative 
approach, keeping t 2 and e 2 /C as the small parameters 
at the same order. Without Ht/c/l/i we are left with 
the ground state of Hbcs Hamiltonian, denoted by | g), 
degenerated by all possible values of <j) and energy E g = 
0 (reference value). All excited states will be denoted 
by |e) with energies E e > 0. Furthermore, we can find 
ground and excited states for Ht included, but now they 
become in general ^-dependent, \gffj with energy E g< f, and 
|e d ) with energy E e( j,. For a given (j) all such states are 
orthonormal and form a complete basis, but states for 
different (j> are not necessarily orthonormal. 


The full eigenproblem reads E\i/j) = H\i/j), where 
the zero order state is | g) with energy 0. Since Ht is 
small, all states \g c / ) ) will also reduce to | g) in zero or¬ 
der. Therefore we can write an ansatz for | iff) in the 
form | ip) = f #(V#)|SV>) + EeV’eWI^)), where the 
latter sum is of higher order. Plugging this form into the 
eigenproblem we get 

E '^\9<t>) + E '^) = ( E g<t> - ^j d l) V’M 

e ' ' 

+ Y1 (' E ^ ~ | £ d l ) V’eM (3) 

with E' = E — Hl/j. Sandwiching it with (g ^| we get 
E 'ip = E g ^ - (d, l> tp){e 2 /C)(g 4l \d l i > \g < p) 

-^{e 2 /2C)( grk \dl\g^) - ^(s^V’ele*)- (4) 

e 

The factor {g^d^g^} is an analogue to the differential 
Berry phase |18j . which can be arbitrarily chosen. For 
our purpose, it is convenient to assume that it is zero. 
Only relative Berry phase for excited states would matter 
but only at high perturbation order. In the lowest order, 
we obtain 

E '^ - E 94> ^ | £$1 - ipie 2 l 2 C){g<t>\dl\g<t>)- ( 5 ) 

The first two terms can be interpreted as adiabatic ap¬ 
proximation m while the last one is a non-adiabatic 
correction to E gc f, = e_ which is the main term we are 
here interested in. The next order term is a bit more 
technical and contains see (T2J. Even higher terms 
get very complicated and very likely already obscured by 
other effects so there is no practical reason to discuss 
them. 

The term {g^ld^lg^) can be evaluated inserting identity 
between derivatives, {g^d^g^) = (g^d^g^,) (g^d^g^,) + 
J 2 e(g<t‘\d<l>\ e <l>)(e 4 >\dci>\g ( t,) since we assume zero differen¬ 
tial ground Berry phase, and from orthogonality be¬ 
tween g and e we get = - E e I(s</>l<^|e0)| 2 = 

— E e \( e <t>\®<t>\9<t>)\ 2 The final form of the equation with 
lowest correction is 

- Eg^ ~ ^ d 2 ^ + V>(e 2 /2C)c(0) (6) 

Here c{<j>) = E e l/e^l 2 , with f e(j> = (e^|9 0 | g^), is the 
(dimensionless) prefactor of the main correction to the 
Josephson energy As we consider the zero temper¬ 
ature case, in the ground state all single-particle states 
with negative energies are fully occupied while positive 
energies are fully empty. The relevant excited states 
(for the terms {g\d\e)) are those with a single excita¬ 
tion, namely one state with positive energy gets fully 
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occupied at the cost of full emptying one of the nega¬ 
tive energy states see Figj2j Note also that for many 
independent channels both terms have to be summed 
over channels, indexed by their transmission t, namely 
Y2 t Egc/, and Y2t c *( < / > )- One can write down c(</>) in terms 
of first-quantized (oc-cupied and em-pty) single particles 
states 


c{4>) = ^2 \{^em\d^\lpoc)\ 2 ■ (7) 

oc,em 

One can also use the adiabatic identity 

teftM = (8) 

because only Ht depends on <j>. This is already sufficient 
to find the lowest contribution in t. In this approximation 
states g and e become ^-independent. Putting E g = 0 
as a reference we are left with excited states represented 
by |d + ,6L) where 6± represent the single-particle state 
in the positive/negative spectrum. Since Ht transfers 
between 1 and 2 the states 6 are in different supercon¬ 
ductors. In momentum space we have k = Asinhd and 
\6±) = (2cosh0) _1 / 2 (e ±6>± / 2 ,±e =F6l± / 2 ) T in the eh basis, 
with e±(6) = iAcoshd and E e — e + (9 + ) — e_(0_). In 
this approximation 


c{<j>) ~ 2 1 2 


dO-dO+ | e (fl + -fl-)/2+i0 _ e (fl_-8 + )/2-i*|2 
(47r) 2 (cosh 9 + + cosh 6L) 2 

(9) 


By introducing variables 2s = 9 + +d_ and 2 w = 9+ — 
we get 


c{4>) -1 2 


dsdw sinh 2 w + sin 2 < j> 
(27t) 2 cosh 2 s cosh 2 w 


( 10 ) 


The integral contains an infinite but phase-independent 
part limited by energy cutoff Cdi V ~ (t/7r) 2 ln(e max /A). 
Interestingly this term appears already in the normal 
case, and can be related to Fermi edge singularity m- 
There are no more any infinities in higher order terms, 
and so this infinity has no physical relevance. The 
phase-dependent part is fortunately finite and equal 
t 2 sin 2 (f>/n 2 . This is our main result. The use of canoni¬ 
cal quantization is valid only if the correction to Joseph- 
son energy is small, i.e. e 2 /2 C -C A. This confirms 
the intuition da that the correction is small at phase 
variation ~ Ec/h ( Eq = e 2 /2 C) smaller than A /h. Re¬ 
member that the quantum operation regime leads also 
to the condition l 2 < 1 or Y2t i 2 ^ 1 in the multi¬ 
channel case. Note also, that in the lowest order only 
A appearing in Josephson energy gets renormalized to 
A' = A + e 2 /2C , 7r 2 ; there are no qualitative differences. 
The qualitative change appears in higher order terms 
~ t 4 but also in the adiabatic ground energy E g ^, so the 
whole description gets anyway complicated J2U]. Note 
that higher order terms will also contain corrections to 



FIG. 3: The exact dependence of c((f>) on <j> £ [0, 7t/2] and t 
referenced at <f> = 0, c is even in (f> and has the period 7r. 


A near the tunneling point (also in E gi p). Nevertheless, 
one can calculate exact c((f>) analytically for all t (not 
only small) assuming constant (bulk) A, which is quite 
lengthy and requires considering all possible excitations 
shown in Fig[2] see details in the Supplemental Material 
Il2l . Such excitations are relevant also in other applica¬ 
tions of ABS m- Here we present this general result 
only graphically in Fig. [3] 

In conclusion, we have found rigorous condition of va¬ 
lidity of canonical quantization of the Josephson effect 
in presence of capacitance, e 2 /C <C A, showing that the 
lowest correction simply renormalizes Josephson energy, 
and constructing a general procedure to find higher cor¬ 
rections. The lowest correction may be useful in fine- 
tuning of quantum energy levels. 
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FIG. 4: The propagation of e-lectrons and /i-oles. We shall use the notation R/L for the inodes at x > / < 0 


Supplemental Material 


A. SYSTEMATIC PERTURBATIVE APPROACH TO JOSEPHSON EFFECT WITH CAPACITANCE 

We shall show how to get at least two lowest order corrections to the adiabatic canonically quantized Hamiltonian, 
outlining the general procedure. We start from ([3]) and Q from the main text. Sandwiching ([3| with (e^, in the 
lowest order, we get 

E e ^e ^ (e 2 /2C')<e 0 |a|'0l^> (A.l) 

Plugging it into the last term of Q we get the desired higher order term. Note also that in the lowest order 
(g<l>\d^ip e \e,p} ~ 2(d0-0 e )(S0l<9</>l e 4-) and g^) ~ 2{d^)(e^\d^\g^). The final form of the equation with two 

lowest corrections is 

E'i> * E g ^ - + V’(e 2 / 2 C) M + ^2 E /eW^eV WiW)) (A- 2 ) 

e e 

with f e<j> = (e 0 1 <9^| g^). 


B. ANALYTIC CALCULATION OF c(0). 
Scattering matrix 


To operate on the continuum states we need to find their scattering matrix. Note that the natural directions of 
electron/hole propagation as in Fig. [4j may be different from in/out-going scattering eigenstates of hscs because the 
direction depends on the sign of group velocity de/dk as depicted in Fig. [5] 

We have 4 incoming and 4 outgoing modes depicted below, indexed by nAa, where n = 1,2 denotes the supercon¬ 
ductor A = R,L (right/left) denotes the side of wave (R = x = 0+, L = x = 0—) and a = i,o denotes i-ncoming and 
o-utgoing mode,see also Fig|U 

The scattering matrix connects the modes amplitudes 


A 


O 



^ A\ Ro \ 


f A\Li\ 

(A Ro \ 

A2 R o 

— Q 

A2L1 

Wo) - 

A-ILo 

— 0 

A i R i 


\A2L0J 


\A21u ) 



= SAi 


(B.l) 


We take 9 > 0. For energies iAcosh# we take 9 Ro = On = ±9 and 9m = 9 l 0 = T&- Then the boundary condition 
in normal scattering matrix ([2]) leads to 

e ^A Ro + e^Am = S e (e s / 2 A Li + e~ 6 ^A Lo ) 
e~°l*A Ro + e e ^A Rl = S h (e~^A Li + e e ^A Lo ) 


(B.2) 
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FIG. 5: The k dependence of direction of BCS eigenstates, based on group velocity dt/dk. 


equivalent to 


so 


We have to solve 


which is equivalent to 


( e*/ 2 ~e~ 9 / 2 S e \ (A Ro \ ( e 9 ' 2 S e -e~ 9 / 2 \ (A Li \ 

\ e -B/2 _ e 0/2 Sh J [ Al J ~ [ e -o/2 Sh _ e e/2 J [ Ar J 


( e 9 ' 2 _ e -B/2 S \- 1 f e e / 2 S e -e~ 9 / 2 \ 
\ e -e /2 - e o/ 2 Sh J [ e -B/ 2 Sh _ e */2 j 


( e 9 ’ 2 -e~ 9 / 2 S e \ (ab\_, 
\e~ 6 /* - e 9 / 2 S h ){c d)- 1 


e 9/2 a - e~ 9/2 S e c = 1, e 9/2 b = e~ 9/2 S e d 
e -B/2 a = e B/ 2 ShC , e -e/ 2 b _ e e/ 2 Shd = J 


Substituting a = e e ShC and b = e e S e d we get ( e 8 Sh — e 9 S e )c = e 9 1 2 and (e 9 S e — e 9 Sh)d = e 9 1 2 
Let us focus on the matrix W = e 9 Sh — e~ 9 S e = ( e 9 Sh — e~ 9 S e ) which appears in 


We have 


( e 9 ' 2 -e- e / 2 S e y l 
\e~ 9 / 2 -e 9 / 2 S h ) 


(e 9 / 2 S h 
V e- 9 ' 2 


-e~ 9l2 S e \ 

-e 9 ! 2 J 


w - 1 


S h - e~ 6 S e = 2 


r sinh 6 it sinh(0 — i<f >) 
it sinh(6* + icf >) r sinh 6 


(B.3) 


(B.4) 


(B.5) 


(B.6) 


(B.7) 


(B.8) 


The determinant is equal 4M, where M = | sinh(0 + iifi )| 2 = sinh 2 9 + t 2 sin 2 cf>. Then W 1 = (AM) 1 (e 9 Sj l — e 9 S']) 
Finally 


4 MS = 



-e~ 9 / 2 S e \ 

—e ®/ 2 J 


(e 9 St-e~ 9 Si) 


( e 9 / 2 S e 
\e~ e / 2 S h 



or explicitly 

/ r sinh 2 (9 — f sinh0sin(</> — iO) — t 2 sin^sin((/> — iO) — rt sinh 6 sin <j) \ 

t sinh 9 sm(<j) + iO) rsinh 2 # rt sinh 9 sin <f) — t 2 sin^sin(0 + i9) 

—t 2 sin (j) sin((f> + i9) — rt sinh 9 sin (j) rsinh 2 0 — t sinh^sin^ + W) 

\ rt sinh 9 sin (f) — t 2 sin<^sin(</> — i9) f sin# sinh (</> — i9) rsinh 2 # / 


(B.9) 


(B.10) 



Andreev bound states 


Apart from the continuum states, there are evanescent subgap solutions of the equation e = h for h = Iibcs + hr 
given by (JTJ) . For evanescent modes e KX i/>, denoting k = Asin/3, |/3| < n/2 we get eigenvalues e± = ±Acos/3 with 
eigenvectors 


ip± 


2 -i/2 


e TM 3/2 \ 

±e^) eh 


Certainly these modes can exist only in the vicinity of the junction point x = 0 forming bound states. The Andreev 
bound states are formed by left (L for x < 0) and right (R for x > 0) evanescent modes. Let us take (3 R = — (3l = f3 < 0 
for the e + and (3 R = —/3l = —ft < 0 for e_. We use the boundary condition at x = 0 to get 


e~ i0/2 A R = e i ^^ 2 S e A L 

e ip ' 2 A R = e~^l 2 S h A L (B.ll) 

for the scattering matrices given by Q. Hence A R = e l ^S e A R = e _z/3 ShAl The condition for existence of solution is 


0 = det(e i/3 S e - e~ ip S h ) = 4det 


f ir sin /3 
t sin(/3 — </>) 


—t sin(/3 + </>)\ 
ir sin f3 J 


(B.12) 


The last determinant is equal — r 2 sin 2 /3 — t 2 sin 2 (3 cos 2 cf> + t 2 cos 2 f3 sin 2 </> which is — sin 2 j3 + t 2 sin 2 (f>. Since it must 
vanish, we get the condition 


sin f3 = ztt sin </> 


(B.13) 


the amplitudes are found by ir sin/Mm = f sin(/3 + </>)Al 2 - Plugging ( B.13| ) we get ±ir svn.(j)ALi = (±t sin </)cos </> + 
cos/3sin</>)Ai2 which lead to ±irA R i = (cos/3±t cos</>)Al 2 - Note that cos f3 = (1 — f 2 sin 2 q i) 1 / 2 = ( r 2 + t 2 cos 2 q i) 1 / 2 > 
\t cos</>| and (cos/3 + 1 cos 4 >)(cos f3 — t cos</>) = cos 2 (3 — f 2 cos 2 </> = r 2 so we have ±?'(cos/3 =F tcoscl)) 1 ^ 2 A L i = (cos/3 ± 
tcos^) 1 / 2 Ai 2 - The normalized solution reads 


A _ o—l/2 ( (1 ±t.COS(f)/COS/3) 1 / 2 \ 

\±*(1 =F tcos</>/cos/3) 1/,2 y 


Rewriting equations (B.ll I in the form 


e~ lf3 A R = S e A L 
e lf) A R = ShA L 


and subtracting the second from the first equation we get 

- 2i sin f3A R = (S e - S h )A L = 2 


0 — tsin< 

f sin </> 0 


and finally, plugging (B.13) and (B.14) we get 


A. R = 


( (1 =F t cost/)/cos/3) 1 / 2 \ 
±i( 1 ± t cos (/>/ cos /3) 1 / 2 / 


(B.14) 


(B.15) 


(B.16) 


(B.17) 


The final wavefunctions have to be normalized by multiplying by \t sin</)| 1//2 . 


Berry phase 

The Berry phase itself will be irrelevant for the lowest orders of perturbative approach. Nevertheless, we shall need 
it because of usefulness in other calculations. We ask what happens when the eigenbasis is parameter-dependent. The 
parameter is </> in our case. For the nondegenerate eigenvalues the differential Berry phase is defined by (^d^ip) = i"f 
where 7 is real because 0 = dq,{ip\ip) = *(7 — 7 *) For degenerate case we collect all the states of the same energy. 
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Assume (it will be anyway our case) that degeneracy is not lifted by changing <p (except some singular points). 
We construct matrix 7 with (ipruld^ipn) = * 7 mn Note that 0 = d<i,(ip m \ip n ) = *( 7 mn — Inm) which shows that 7 is 
Hermitian. We can find individual Berry phases by diagonalizing 7 . Then eigenvalues will be differential Berry phases. 

Let us calculate first the Berry phase for Andreev bound states. The derivative of spatial part is real so it cannot 
give any contribution to Berry phase. Acting on e ±l/3 / 2 the net result is also zero as both signs appear with equal 
weight. Finally acting on A R / L we get real numbers, so again they cannot give any Berry phase. The conclusion is 
that the phase convention we used in construction of Andreev bound states representatives gives zero Berry phase. 

For scattering states we have the total state combined from incoming wave and scattered outgoing wave tp a = 
Aa + E c il>ocSca where a,c = R/L , 1/2. Only S can depend on <p so the matrix 7 is given by 

*7 ba = X! I ^o^dfSca (B.18) 

The overlap of Ri/Ro and Li/Lo is ~ i/sinh0, so it does not contain any <5 term (the i/o states have opposite 
momentum) so the term vanishes in the continuum limit (infinite length of the junction). Therefore i'y = S'd^S/2 
because the outgoing waves are halves of full waves. The final calculation gives 2 M 7 = 


t sinh 9 


/ —t cosh 9 ir cosh(0 + icp) 
—ir cosh(0 — icj>) t cosh 9 

0 i cos <j) 

\ —i cos <p 0 

We shall also need *7 = Sd^S^/2 which is 2M 7 = 

/ — t cosh 9 — ir cosh (0 + icp) 

ir cosh{9 — icp) t cosh 9 

0 —i cos 6 


t sinh 9 


V 


l cos ( 


0 


0 

—i cos cp 
t cosh 9 
—ir cosh (0 + 


0 

i cos <p 
t cosh 9 


i<t>) 


% cos cp 

0 

ir cosh (0 — icp) 
—tcosh 9 


—* cos cp 

0 

—ir cosh (0 — i(j>) 


ir cosli(0 + i(j>) —t cosh 9 j 


(B.19) 


(B.20) 


Contribution to c((j>) from ABS 


We shall calculate all contributions to c defined by ([7|) by inspection of all relevant excited states, depicted in Fig{2] 
Let us first calculate A = where ip± is one of Andreev bound states, while ip is any eigenstate. If ip is the 

opposite Andreev bound state there is no contribution from spatial part |«| 1 ' 2 e~l rex l for the same reason as the Berry 
phase vanishes. If the derivative act on e/h components then the contribution cancels from orthogonality of Ai/ 2 - 
The only nonzero term is when acting on A 1 / 2 . The final result is 

{ip±\d ( j,\ipzf.) = ±it 2 r sin 2 0 / 2(1 — f 2 sin 2 <p) (B. 21 ) 

and ca = t 4 r 2 sin 4 </>/4(1 — t 2 sin 2 <p) 2 

The relevant cross terms between ABS and scattering states are when their energies have opposite signs. Therefore 
the sign ± in ABS will determine the energy (its sign is the sign of —j3 = =Ffsin/>). We are interested in the vector 
(ip±\(\ipiRo), \ip2Ro), | ipiLo), \' l p2Lo))l because act only on scattering matrix, d^S, accompanied by outgoing states 
and we can insert at the end (since the relevant quantity is Hermitian square). The result is 


— 2i\t sin^l 1 / 2 sinh ((0 + if))/ 2 ) 

( 2 cosh 0 ) 1 / 2 (|tsin</| — *sgn/3sinh0) 

((1 =F tcoscp/ cos /3) 1 / 2 , =F *(1 ± tcoscp/ cos /3) 1//2 , — (1 ± tcoscp/ cos/)) 1 ^ 2 , H( 1 =F tcoscp/ cosf)) 1 ^ 2 )^ 
The final element cbb = Y/ a \(M^ci>\' l Pa)\ 2 is equal 

1 1 3 sin cp\ sinh 2 9( cosh 9 — cos [3 + 2 r 2 sin 2 cp/ cos f3) 

(cosh 9 — cos /3)(cosh 9 + cos /3 ) 3 cosh 9 

Its total relevant contribution is cb = / 0 °° cbo cosh 9d9/2n which is 


|t 3 sin</| /2 ^cos 2 /3 
87tcos 4 /3 V sin/3 


2 cos 3 /3 — 7 r| sin /3\ cos 2 /3 + cos 2 <f> ( 7 r| sin^| + 


2/3 cos 2/3 — sin 2/3 
sin 3 /3 


For small t (/3 —>■ 0) the leading terms are |f 3 sin </>|(1 + 2 sin 2 </>)/ 67 t — t 4 sin 4 (p/8 


(B.22) 


(B.23) 


(B.24) 
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Contribution to c(</>) from scattering states 


We are interested in terms cg d = J2ab I(^o,|<9^|i9 &)| 2 where 1 9 a ) is a scattering state with energy — Acoshf? while |t?b) 
is a scattering state with energy +A coshth By analogous considerations as with ABS, the derivative effectively gives 
(d ( f,S)S^ = —*27 Therefore c is equal 4Try 2 (9)Q where Q a b = X) c ($ O a|0c)(0c|$oi>) The middle part is 


£i<w«y = £ \0ic){e ic \ + \e od )(ss^) dc (9 oc \ + \0oc)s cd (e id \ + h. c . 

c c dc cd 

From unitarity of S we get 

£iWci= £ \0xc)(9xc\ + \9oc)S cd (0id\ + h.C. 


(B.25) 


(B.26) 


C’,X=l,0 


cd 


Therefore Q = BiB\ + B„Bl + BiS^B\ + B 0 SB\ where B x = (9 0 \9 X ) with x = i,o. We find 


B, = 


isinh((i? — 9)/2) /0 / 

(cosh 6 cosh 9) 1 / 2 (sinh 9 — sinhd) I 0 


, B 0 = 


—i sinh((i? + 9)/2) // 0 

(cosh 9 cosh i ?) 1 / 2 (sinh 9 + sinh 9) \0 —I 


where / is 2 x 2 identity matrix. Plugging in all matrices we get 


cg d = 


(2t sinh 0 sinh d ) 2 (cosh f? cosh + 1 — 2 cos 2 (f> — t 2 sin 2 <j>) 
cosh 9 cosh 9 (cosh 9 + cosh 9) 2 (sinh 2 9 + t 2 sin 2 0)(sinh 2 9 + t 2 sin 2 </>) 

2 


(B.27) 


(B.28) 


The final relevant integral is cc = / 0 °° d sinh 9d sinh 9cg d / (2tt) 

We can divide the integral into divergent and convergent parts c c = c d i V + c con with 


Cdiv 


dOdd 4t 2 (cosh 9 cosh d + 1 — 2 cos 2 0 — t 2 sin 2 0 ) 


(2tt) 2 


(coshf? + coshd ) 2 


and 


d0d9 4 1 4 sin 2 </>(cosh 9 cosh 9 + 1 — 2 cos 2 <j> — t 2 sin 2 <f>) (t 2 sin 2 <fi + sinh 2 9 + sinh 2 9) 
(27r) 2 (cosh# + cosh 9 ) 2 (sinh 2 9 + t 2 sin 2 <^)(sinh 2 9 + t 2 sin 2 cj>) 


(B.29) 


(B.30) 


In the first integral we make substitution 2 A = 9 + 9, 2 B = 9 — 9. The calculation needs a cutoff in the first term at 
.4 jff a x and then gives 7 r c d iy — t ^^uiax t cos (j) t sin cJ)J 2 Aow Ccon — C-conl + Cconl with 


Cconl — 


Ccon2 — 


f°° 2 1 4 sin (j)d0dd cosh 9 cosh d + 1 — 2 cos 2 (j) — t 2 sin 0 

Jo 7 r 2 (sinh 2 9 + t 2 sin 2 4>) (cosh 9 + cosh 9 ) 2 

30 2 1 6 sin 4 cj>d9d9 {cosh 9 cosh 9 + 1 — 2 cos 2 (f> — t 2 sin 2 <j>) / 7 r 2 


/o (cosh 9 + cosh 9) 2 (sinh 2 9 + t 2 sin 2 cj>) (sinh 2 9 + t 2 sin 2 <j>) 
We shall integrate over 9 with help of identities 


/»oo 

/ d#(coslr# + coshd)^ 1 = 9/ sinh 9 

J o 

f°° ,,, . . , q ._2 9 cosh 9 - sinh 9 

/ d9(cosh9 + coshi?) =-=- 

Jo sinh 9 

f°° ... . . _ .. o 9(3 + 2 sinh 2 9) — 3 sinh $ cosh d 

/ d9 (cosh 9 + cosh 9)-* = v 

Jo 


(B.31) 


(B.32) 


2 sinh 5 9 


the integration gives 


Cconl — 


2 ddt 4 sin 2 0 / 7 T 2 f d cosh d 


sinh 2 d + t 2 sin 2 0 \ sinh d 


— (sinh d + t 2 sin 0 + 2 cos 2 0 ) 


, 9 cosh 9 — sinh 9 \ 
sinh 3 9 J 


(B.33) 
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We shall integrate C c 2 over 0, noticing residua in i/3, i( 27 r — /3), i( 7 r ± /3), and ±d + in. 

2t 6 sin 4 </>dd 


Ccon2 — 


7 r 2 (sinh 2 d + sin 2 /3) 

(n — /3) (cos/3 cosh d + cos 2 /3 — 2 cos 2 <A) /3(cos 2 /3 — cos /3 cosh d — 2 cos 2 </>) 


2sin/3cos/3(coshd + cos/3 ) 2 


2 sin /3 cos /3(cosh d — cos /3 ) 2 


1 — dcothd 2 cos 2 </>(! — dcothd) — d sinhd coshd 4 cos 2 </>d sinh 1 ? cosh d 


sinh d 


+ 


sinh 2 d(sinh 2 d + sin 2 /3) 


sinh 2 d(sinh 2 d + sin 2 /3 ) 2 


The integration over d needs noting that sinh 2 d + sin 2 /3 = (coshd + cos /3) (coshd — cos /3) and further ir 
simplicity we assume (3 £ [ 0 , 7 r]) 


nOO 

Pl(/ 3) = / dd(coshd + cos/3 ) -1 = /3/sin/3 

Jo 


and by derivatives (sin 1 (3d/df3) 


P 2 (/ 3 ) = / dd(coshd + cos/3) = 


_ 2 sin /3 — /3 cos /3 


/»oo 

Pz{/3) = / dd(coshd + cos/3 ) -3 = 

Jo 


sin 3 /3 

_ 3 /3(3 — 2 sin 2 /3) — 3 cos/3 sin/3 


2 sin 5 /3 


Also 


Pl(/ 3) = / dd(sinh 2 d + sin 2 /3) 1 = (Pi(7r — /3) — Pi(/3))/2cos/3 = (7r/2 — /3)/sin/3cos/3 
3o 

/»oo 

A(/3) = / d^(sinh 2 $ + sin 2 /3) -2 = —P{(/3)/2sin/3cos/3 

Jo 

/>oo 

PoiP) = / dd(sinh 2 d + sin 2 /3) -2 = — P 2 (/3)/4 sin/3 cos/3 
3o 


A more complicated integral 

QiOS) = 

And its derivatives 

Q 2 (/3) = 

W) = 


d cosh ddd 


/o sinh d(cosh d + cos/3) 


2/3 2 cos /3 + 7 T 2 (1 — cos /3) 
4 sin 2 /3 


Z 100 dcoshddd 7 r 2 (2 — sin 2 /3 — 2 cos /?) + 4/3 sin/3 cos/3 + /3 2 (2 sin 2 /3 — 

0 sinh d (cosh d + cos/3 ) 2 4sin 4 /3 

Z 100 dcoshddd 


0 sinh d(cosh d + cos/3 ) 3 


7 r 2 (1 cos/3) + 2 cos /3 sin /3 + 6/3sin/3(sin 2 /3 — 2) + /3 2 (2 cos/3 — sin /3 cos /3 + 2 cos/3 sin 2 /3) 


4 sin 6 /3 


/»oo 

/ dd(sinh 2 d + sin 2 /3 ) _1 = 4 
do si 


tt/2 - /3 


sin /3 cos /3 


Another important integral is 


*,(/>) = r«- 2 = / 

do sinh d(sinh 2 d + sin" /3) do 


00 dsinhdcoshd. . , _ 2 . . . , 2 „ . 2 a\-ii 

dd- : -[smh " d — (sinh" d + sir /3) ] 


sin /3 


Differentiating Pi sin" /3 over /3 we get 


r °° 2 sin/3 cos/3 sinh dd cosh d 

ad- ; - 33 ——- 

- -„2 .a , „:„2 o \2 


/0 


(sinh 2 d + sin 2 /3) 


(B.34) 

itegrals (for 

(B.35) 

(B.36) 

(B.37) 

(B.38) 

4) 

(B.39) 

(B.40) 

(B.41) 



12 


Integrating it by parts we get 


d&- 


sin /3 cos P 


= (tt/2-/3) 


Jo (sinh 2 1 } + sin 2 /3) 2 
from the previous integral. We note also that R\{P) sin 2 p —>• 0 for p —>• 0 so 

R\{P) = p{n — /3)/2 sin 2 p = (Qi (tt — /3) - Qi(/3))/2cos/3 

Its derivatives 

$ cosh $ 


MP) = [ 

Jo 

*3 09 ) = / 

Jo 

It gives c co „i7r 2 /t 2 equal 

sin 2 /3 + /3 2 — 7 t|/ 3| + 2 cos 2 (f> 
and c con2 7r 2 /2f 6 sin 4 <f> equal 


dti- 

c 

dtf- 


= — R[ {pi)/2 sin /3 cos p 


0 sinh $(sinh 2 •& + sin 2 P) 2 

^ C0Sh ^ 2 ^ = ~ R 2 0 3 )/ 4 sin P cos P 


0 sinh i?(sinh 2 1 ? + sin 2 /?) 
7r| sin p — P cos P\ 


sin 2 P cos P 


2 cos P — 2/3/ sin /? + 2 


P(pcosp — sin/3) 
sin 2 P cos P 


with a equal 


a(/3) + 6(/3) cos 2 0 

Pi(tt - /3) /Pr(7r - p) P 2 (/3) , Pr(/3) 


+ - 


4 cos 2 p 2 cos /3 4 cos 2 p 

Pi(p) ( PliP) P2{*-P) Pi(tt — /3) 

4 cos 2 P 2 cos P 4 cos 2 p 


~MP)- 


Ri{P) - Pi(P) — 1/2 


sin 2 P 


and b equal 

Pr (tt - /3) / A (tt - /3) P 3 (/3) P 2 (/3) P, (/3) 


COS P 

PliP) 


8 cos 3 P 2 cos P 4 cos 2 P 8 cos 3 P 
p i{P) , Po{n-p) P 2 (tt-/3) + Pi(tt-/3) 


(B.42) 


(B.43) 


(B.44) 


(B.45) 


(B.46) 


cos /3 \ 8 cos 3 /? 2 cos /? 4 cos 2 /? 8 cos 3 (3 


-4 P 3 (/?) + 


(B.47) 

2P li 5 — 2Pi(/3) — 1 , 2P 2 (/3) — 2P 2 (P) 


sin 4 P 


sin 2 P 


In the lowest order c con ~ —|< 3 sin^|(l + 2 sin 2 4>)/3tt + 2t 4 sin 4 0/7T 2 +f 4 sin 4 0/4 so the f 3 terms from ABS cancel! 



